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ABSTRACT 


The aim of this article is to determine the stability characteristics of a Rayleigh 
layer, which is known to occur when the fluid above a flat plate has a velocity im- 
parted to it (parallel to the plate). This situation is intrinsically unsteady, however 
as a first approximation we consider the instantaneous stability of the flow. The 
Orr-Sommerfeld equation is found to govern fixed downstream wavelength linear per- 
turbations to the basic flow profile. By the solution of this equation we can determine 
the Reynolds numbers at which the flow is neutrally stable; this quasi-steady approach 
is only formally applicable for infinite Reynolds numbers. We shall consider the large 
Reynolds number limit of the original problem and use a three deck mentality to de- 
termine the form of the modes. The results of the two calculations are compared, and 
the linear large Reynolds number analysis is extended to consider the effect of weak 
nonlinearity in order to determine whether the system is sub or super critical. 


* Research was supported by the National Aeronautics and Space Administration under NASA contract 
No. NAS1-19480 while the author was in residence at the Institute for Computer Applications in Science 
and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23681 
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§1 Introduction 

The main thrust of boundary layer stability calculations have involved spatially 
evolving problems. Relatively little work has gone into the stability of temporally 
developing layers, one such situation is the topic of this paper, that of an impulsive 
flow over a flat plate. 

The calculation to determine the linear stability of the system considered here is 
initially attempted using a quasi-steady approach where the flow is frozen at a certain 
time. This yields an Orr-Sommerfeld equation, the solutions of which were discussed 
for a variety of flows by Tollmien (1929) and Schlichting (1933), using homogeneous 
assumptions in the temporal and downstream coordinates. Unfortunately the true 
system governing what we shall refer to as Tollmien-Schlichting waves is not uniform 
for a general boundary layer. An attempt to remedy this deficiency was proposed for 
the Taylor-Gortler problem by Smith (1955). These ideas were extended by Bouthier 
(1973) and Gaster (1974), where extra terms were added to appreciate for the boundary 
layer’s growth; although it was noted by Gaster that this only amounted to a successive 
approximation technique. It was in the work of Smith (1979 a) that a more formal 
asymptotic treatment of the problem was proffered for large Reynolds numbers. Smith 
used a triple deck structure to describe the linear development of fixed frequency 
Tollmien-Schlichting waves within a growing boundary layer, specifically a Blasius 
layer. It was found that the boundary layer growth altered the higher order terms in 
the expansions of various wave characteristics within this structure. The inclusion of 
these terms may be stabilising or destabilising depending on which quantity is taken 
to be representative of the perturbations intensity. 

Some studies have been made of the stability characteristics of unsteady bound- 
ary layers, although most of the cases have involved temporal periodicity. In these 
problems either a quasi-steady approach or Floquet theory can be used. In Kerczek 
& Davis (1974) the linear stability of a Stokes layer on a flat plate is discussed, in this 
problem, waves that develop quickly compared to the temporal scale of the basic state 
are considered, as is the case in Cowley (1986). In Hall (1975) a study was made of 
unsteady flows about cylinders, and in Seminara k Hall (1976) the Taylor problem is 
considered where the inner cylinder is taken to oscillate periodically with time. The 
work of Di Prima k Stuart (1975) regards the Taylor problem for eccentric cylinders 
(which corresponds to a journal bearing). In Otto (1993) the stability of a boundary 
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layer at an infinite cylinder is discussed. The cylinder starts to spin up and vortices 
are generated by an axisymmetric. element of wall roughness; these may lead to some 
kind of transition. In Duck (1990) the unsteady triple deck equations were solved for 
the flow over a small hump on a flat plate. 

In the current problem the basic state depends on the normal and temporal coor- 
dinates, so it is natural to solve an unsteady problem. This situation is fully parallel 
so there is no need to make any assumptions about the spatial evolution of the basic 
state. After a time to? we may define the boundary layer Reynolds number as, 

„ _ Uoo{vuy _ UooIq 

•/teO i j 

V 1/7 

based on the viscous layer’s thickness. The velocity imparted to the fluid is Uoo } 

and its kinematic viscosity is v . In the viscous layer (which may be shown to have 
1 

thickness (vt) 1 ), the basic state is found to be similar in time and is given by an error 
function. This is accreditable to the large Reynolds number assumption; at finite 
Reynolds numbers the unsteady equations need to be solved. We wish to consider the 
susceptibility of this profile to Orr-Sommerfeld modes that evolve on short temporal 
and spatial scales. As mentioned previously the approximation is made that we shall 
‘freeze’ the basic profile. In an infinite Reynolds number flow this is formally correct 
as there is no way information can propagate; but in the Orr-Sommerfeld method 
the Reynolds number is retained as a finite parameter. In this paper we do not 
wish to discuss the legitimacy of the Orr-Sommerfeld approach, it is just used as a 
first approximation to the solution of the linear stability problem at finite Reynolds 
number. For each time to beyond a certain point (this corresponds to Reynolds number 
R e o), this problem may be solved to obtain a real wavenumber and phase speed such 
that the wave is neutral. Above this critical Reynolds number there are two neutral 
values, and as the Reynolds number increases still further the behaviour of the wave 
having the smaller wavelength is governed by a five deck structure, the critical layer 
becomes distinct from the wall layer, due to the size of the phase speed. Our interest 
is confined to the lower branch where the structure may be described using three 
decks. An ample description of this structure can be found in Stewartson (1974). The 
bulk of the boundary layer is referred to as the main deck, which at leading order is 
inviscid and linear. In the proximity of the wall there is a fully viscous layer, this 
region contains the critical layer of classical Orr-Sommerfeld theory. Finally there 
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is an upper deck which is external to the conventional boundary layer and enables 
matching with the outer potential flow solution. A multiple scales technique is used 
to resolve the growth of the boundary layer with time. 

The temporal growth of the layer produces corrections to the eigenproblems for 
the wavenumbers and phase speeds at a higher order. The triple deck solutions may 
be compared with the Orr-Sommerfeld calculation in the regime of common validity, 
that is where i? e o >> 1. In fact the corrections caused by the growth of the layer to 
the ‘steady’ triple deck scenario occur in the third order equations, one order earlier 
than those in Smith (1979 a). This is due to the different temporal and spatial scales 
over which Tollmien-Schlichting waves develop. 

The results of this calculation may be viewed in another sense. For a perturbation 
wavenumber a, at each time to we have a frequency fl. When the imaginary part of 
becomes zero, the wave is instantaneously neutral. In the neighbourhood of this neutral 
time we may develop a weakly nonlinear theory. A Stuaxt-Watson amplitude equation 
can be derived, and using this we may determine whether the perturbations are sub 
or super-critical, that is whether the nonlinearity has a destabilising or stabilising 
effect on the situation. We shall also consider the interplay between the unsteady and 
nonlinear effects, as was discussed for the Blasius problem by Smith (1979 b) and Hall 
& Smith (1984). 

The procedure adopted in the remainder of this paper is as follows, in section 2 
the basic flow is derived and the quasi-steady problem is formulated. In section 3 the 
linear triple deck problem is formulated and solved. In section 4 the weakly nonlinear 
problem is derived and the nonlinear amplitude equation is given. Section 5 contains 
details of the numerical solution of the linear problem. Section 6 includes discussion 
of the effect of considering certain small but finite disturbances, and in section 7 some 
brief conclusions are drawn. 


§2 Derivation of the basic flow and the formulation of the quasi— steady problem 


The problem is non-dimensionalised in the usual way and the unsteady three di- 
mensional Navier-Stokes equations in cartesian coordinates become (with the temporal 
scaling t ~ L/Uoo for the Orr-Sommerfeld modes, where L is a typical lengthscale, and 
t ^ L 2 /v for the basic state), 


du du du du dp 1 , 

m +u Tx+' , d^ + w Tz = -5i + T7 u ' 
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dv dv dv dv dp 1 _ 2 
Ik +U fc +V dy +W fc = ~d^ + W e W V ' 

dw dw dw dw dp 1 2 

^ + u a^ + v d^ + w d^ = -d-z + W t v w ' 

du dv dw 
dx ^ dy + dz ’ 

where R e is the Reynolds number defined as U^L/v, and V 2 is the Laplacian operator 
in cartesian coordinates. The choice of the lengthscale L is discussed in section 5, but 
it shall be taken to be the displacement thickness of the layer, as in Jordinson (1970). 
In this section disturbances shall be considered to be two dimensional to facilitate 
comparison between this work and that of Smith (1979 a), to this end it is assumed 
that djdz = 0 and w = 0. 

We require that the unsteady terms and viscous terms balance in the viscous layer, 

x 

this implies that this layer is of thickness (vt )* . In the non-dimensional coordinates 
this corresponds to . Thus we introduce a boundary layer similarity variable ?/ given 

by 

’"ST - (!U) 

There is no flow in the normal direction, as can be seen from the equation of continuity 
at leading order. Hence the basic state is governed by 

d 2 u B du B 

+ 2r/^^- = 0, 


dif 

along with the boundary conditions 


dv 


u B =0 


at 


V = 


and 


u B 


as 


CO. 


The solution of his system is given by 





o 


this is an error function. 

Now in the context of an Orr-Sommerfeld equation we consider the instant sta- 
bility of this profile, hence (2.1) now becomes r/ = y/2tQ, where t 0 is a constant. We 
suppose that any perturbation to the flow is independent of the order one temporal 
scale. 
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The perturbations are considered to be travelling waves which evolve on the scale 
L, thus the disturbance is proportional to exp(t^(x — ct)). The disturbance is taken 
to be infinitesimal so that the resulting analysis is linear. As the perturbation is two 
dimensional, a streamfunction formulation can be exploited, such that 


U = 


dtp 

V 


v = - 


c )tp 
dx 


It should be noted that terms proportional dtp /dt — 0 have been ignored, so that the 
equation may be solved as an ordinary differential equation. The governing equation 
is the conventional Orr-Sommerfeld equation namely, 


(ub - c ) 



_ Pub = 1 

dy 2 i8R e 





tp. 


( 2 , 2 ) 


The boundary conditions imposed are those of no-slip at the wall, and the decay of 
the disturbance as y — > oo, this implies the wave is confined to the boundary layer. 
This eigenvalue problem can be solved numerically employing the method discussed in 
Malik, Chuang,& Hussaini (1982). The method uses a two-point fourth order compact 
finite difference scheme based on a Euler-Maclaurin method. A stretched grid was used 
so that resolution could be retained without using a prohibitive number of points. The 
calculations were checked by finding the eigenvalues of the adjoint system to (2.2). The 
critical values were found to be 


R ec = 2968.4, n c = 0.13340, 6 C = 0.4255, c = 0.3135, 

where = 8c. This point represents the smallest value of the boundary layer Reynolds 
number at which the flow situation is linearly unstable to travelling waves of the type 
considered herein. Plots of 8 and c against R e are given in figures 1 and 2 respectively. 
These calculation were affected for <o = 1/4, so that y = r) at that instance. 
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Figure l:Streamwise wavenumber versus R e , Orr-Sommerfeld calculation 



Figure 2:Phase speed versus R e , Orr-Sommerfeld calculation 

§3 Formulation and solution of the triple deck problem in the linear regime 

In this section the structure of the the linear disturbances is described using triple 
deck theory. As mentioned previously R t is assumed to be large; and e = Rj 5 <C 1 is 
introduced. In classical Orr-Sommerfeld theory the typical wavelength of the neutrally 

I 

stable modes increases proportionally to J?| as R e 1, see Stuart (1963). Stuart 
confirms this asymptotic property with a Reynolds number based on boundary layer 
thickness. Three distinct regions are considered, namely the upper, main and lower 
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decks, scaled normally on e 3 , e 4 , and e 5 respectively. Further temporal and spatial 
variables are introduced, that is T and X, given by, T = e~ 2 t and X = e~ 3 x. The 
solutions are considered in the three decks in terms of expansions involving powers of e 
and In e, as in Stewartson (1974). The disturbance is now considered to be proportional 
to 


E — exp 



The disturbance is taken to have slowly varying frequency and constant wavenumber, 
so that 0 and Q are given by 


0 — 0\ + e $2 + e 2 In e 03 £ + c 2 #3 + • • • , (3.1a) 


and 

= fii(t) + 5 ^ 2(0 T In l(^) T c 2 f^3 (0 T ■ ■ ■ • (3.16) 

A multiple scales form for the temporal derivative is introduced, 

d_ d_ _ 2 d_ 
dt~* dt +e dT' 

Notice, that to determine the effects of unsteadiness the expansions need only go as 
far as the fourth order terms not the fifth order terms as in Smith (1979 a), this occurs 
due to the difference in scalings associated with the temporal and spatial variables. 

§3.1 The Main Deck 

The transverse coordinate in this layer is Y = 0(1), where Y = e~ 4 y. The main 

_ 1 

deck corresponds to the conventional R e 2 boundary layer. The basic flow in the layer 
is given by the error function as shown in section 2. It should be noted that there is 
no normal component of this velocity, as can be seen by considering continuity in this 
layer in the absence of any disturbance. It may be shown that the basic flow has the 
properties that, 


u B {t,Y) ^ \(t)Y + \ 3 (t)Y 3 +0(Y 5 ) as Y -* 0, (3.3a) 


u B (t,Y) — > 1 as Y — > oo, 


(3.36) 
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where it is known that X(t) - Xt 2 , with A = l/y/n. The perturbation to the basic 
state in this layer takes the form, 


« = [«l + eu 2 + e 2 In eu 3L + e 2 u 3 + . . .] £, 
v = [fuj + e 2 v 2 + ( 3 In eu 3 L + e 3 u 3 + ...] E, 

P = [ e Pi "b f 2 p2 + £ 3 In tpzL, + f 3 p3 + • • •] E. 

Notice that the logarithmic terms (characterised by subscript L), occur at one order 
higher than in the corresponding expansions of Smith (1979 a). This is due to the 
form of u B as Y — * 0 given by (3.3 a). This form of the disturbance is now substituted 

into the governing equations, and nonlinear terms are neglected. The operators M\ 

and M 2 are introduced, and are given by 

Mi (u,v) = iO\u + M 2 (u,v) = iO A u B u + 

where these are identical to those used in Smith (1979 a). 

The equations in this layer are conceptually identical to those of Smith (1979 a), 
except the order at which the logarithmic terms occur. In this deck there is no effect 
from the multiple scales form to the order considered, so it is sufficient just to quote 
the equations governing the logarithmic terms, 


Ml (u 3 l,v 3 l) + i&3LU\ — 0, M 2 (u 3 l,W3l) +UBi0 3 LU\ = 0, 


0 = 


dpi L 

dY ' 


The boundary conditions for the disturbances in this layer are given by matching with 
the lower and upper decks, the solutions for this layer are considered in section 3.4. 

§3.2 The Lower Deck 

In this deck the order one variable y is introduced, defined by y = e~ 5 y. The lower 
deck contains the wall layer and critical layer of classical Orr-Sommerfeld theory. The 
wave in this deck expands as 


u = [U] + £ U 2 + e 2 In tU 3 i + e 2 £/ 3 + . . .] E, 
v — [f 2 ki + ( 3 V 2 + e 4 In eV 3 L + £ 4 V 3 + . . .] E, 

P = [f-P] + f 2 P2 + e 3 In eP 3 i + e 3 P 3 + . . .] E. 

The basic flow is given by u B ~ eyX(t) + e 3 y 3 A 3 (t) + ••• as can be seen from the 
form of (3.3 a). From the normal momentum equation it may be shown that the 
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Pj' s O' — 1,2,3) are independent of y, but Uj and Vj are functions of y and t. Again 
following Smith (1979 a) the operators L\ and L 2 are introduced, and are given by 


L i (u,v) = i9-[U + 


dv 

dy 


L 2 (u,v,p ) = i(\yO] - + \v + i9ip - 

dy 1 


Again it is sufficient only to quote those equations that differ from those of Smith 
(1979 a), 

L\ (U 3 L, V 3 1 ) 4- iQziUx = 0, 


£2 {U3L1 13L, Pzl) — + \yi 0 3 iU\ + iOziPx = 0 , 

■^2 (U 3 , V 3 , P 3 ) — i (Q 2 U 2 + £l 3 U\) 4 Xyi { 62 U 2 + QzU \ ) + i ( 0 2 P 2 4- O^Px ) 

+ iXsSxUxy 3 + 3y 2 A 3 VI 4 — — - = 0. 

at 

The boundary conditions in this layer are those of no-slip at the wall and matching 
with the main deck. 

§3.3 The Upper deck 

In this third deck a coordinate y is introduced, defined by y = e 3 y, where y is an 
order one quantity, and the disturbance takes the form 


u — + c* u 2 + e 3 In €U 3 l ~b c 3 M 3 + . . .] E , 

v = [ev] + e 2 v 2 + e 3 In ev 3 L + ^v 3 + ...]#, 

P = [tpi + <?Pi + e 3 In epzi + e 3 p 3 + ...]E. 
Here the basic velocity is almost the uniform flow solution that is 


ub = 1 + 0 ( e 3 ). 

In this layer all the flow quantities are functions of the slow time scale t and y. The 
operator D is introduced, defined by 

After some algebraic manipulation to remove the velocity components the governing 
equations for the pressure are obtained as, 

D(Pi) = Z?(p 2 ) = 20i d 2 pi , 

E ( p 3 L ) = 2#i03LPi , D (^3) = 26] (d 3 p] 4- 62P2) + 6 \p\ . 
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The boundary conditions are those of matching with the main deck as y — > 0 and 
the requirement that the solution remains bounded as y — > oo. Having set up these 
equations they are now solved in each layer and the eigen problems are derived. It is 
worth stressing at this point that there is no effect in the upper deck equations from 
the unsteadiness (to the order considered). The unsteady derivatives are of a higher 
order unlike the spatial derivative which occur in the diffusive terms, from 





§3.4 Solutions of the main deck equations 

The solutions of the first order equations in the main deck can be shown to be 


v\ = —i6\ A\ub-, 


A ® UB 

U| 


P\ = P \ , 


where P\ = P\ (t) and A\ = A\ (t) are unknown slowly varying amplitude functions. 
The pressure in this deck will be matched with the pressure in the lower deck. Solutions 
here are the same as those given in Smith (1979 a). The integrals H defined in Smith 
(1979 a) are functions of x and F, whereas here they are functions of t and F, also 
obviously ub is the error function, rather than the solution to the Blasius equation. 
The logarithmic solutions are given by, 


— f A ^ 1 ^ 3 l \ dug 

U 3L - 3L ^ J dy , 

V3L = —iOi Azluq, Pzl = Pzl> 


§3.5 Solutions in the lower deck 

To determine the behaviour of U\ in this deck, we differentiate the leading order 
equation of continuity with respect to the decks transverse coordinate, that is y, and 
eliminate V \ . It is required that the disturbance is bounded as y — > oo, thus the 
solution for is found to be given by 

V h = Bi Ai(f), 

where 
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Ai(£) is the Airy function and B\ is an unknown function of t. Applying the no-slip 
condition at y = 0 in the leading order equation, it is found that 


Bi Ai(,A 3 = i9\P\, 


(3.4) 


2 

where subscript 0 corresponds to evaluation at £ = £o = -dli , and a prime 
denotes differentiation with respect to £. The matching condition on the leading order 
disturbance as y — » oo (therefore as |£| — ► oo), is given by U\ — ► AAi, thus it is found 
that 

B\k = \A \ , 


where 

OO 

K = J Ai («) dl 

£<> 

This yields the first relation between the pressure and the negative displacement, Pj 
and A\ respectively, obtained by eliminating Bi from the equation (3.4). Note that 
these are effectively the classical critical layer equations, since, as noted previously on 
the lower branch of the stability curve the critical layer coincides with the wall layer. 
As can be seen from Smith (1979 a) similar techniques may be used to solve for the 
higher order equations. The same solutions as Smith (1979 a) are obtained except the 
unknown amplitudes are now functions of t, rather than x. The equation to be solved 
at the next order is 


d 3 U 2 t dU 2 

d? * di 


—iQ, 2 A' 


>dl h 

' 3 


dt 


+ i A ( £A 3 



A 3 02 


dU^ 

~df' 


which yields the solution 

U 2 £ = B 2 Ai —0 2 d^ ^ B\ A\" Ai^A 5 . 

O 


The function U$l satisfies 

&u 3L r au 3L 

^5--f-sT = - ,n “ A 5 

which gives the solution, 


dU ' j. \f 

W \ 


fA-4 + 


A0j 


A~h, L dUi 


di ’ 


P 3 Lf = B 3L A\ + ^SsLd-'Bi Ai m + *A-3 PiAi'A’i 
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At next order the function U$ satisfies 


d 3 u 3 dU 3 . OA -ia U 2 . 0 A -*dUi 

(~zr=-i\l 2 A 3 -r--zfi 3 A 3 


d( 3 


d( 


d( 


**-' (<*“*&) *^) 


i 9 
+ A 3 — 


a / a 


t_{ AHh fttu AA d_' 
i 3 0\ A y fii A J 


U i 


+ -^-A-ifVi + 6 ^.* 


iA^i 

3 


(w.r 


+ £e< 


i=o 

which yields the solution 


dUi 

d( ’ 


dU : 


+ */«), 


i=0 


where, 


5Ai <ii d 3 Ai 

/(0 -““a^ + y afr + a i 
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<9 7 Ai 3^Ai „5Ai 
+ 2 


de 2 ^ ■ - di 

+ ui A ~ 5 (w ~J J Aid(di ) ~ 7 Hy/ AWe ’ 

where the limits of the integrations are £ to oo. In this solution the terms included in 
the j3 summation include those of Smith (1979 a) and the unsteady terms. The terms 
included in the function /(£) are those associated with A 3 , that is the second order 
mean flow quantity. 

The /3’s are given by, 

A> = #3, A = oui-02 + <*12-^1 — -r-B\ac\\ + A 

"l A 


. A » dB\ 2 A t _ 

A_- 2 ~“ar +Q " B ' ~7 ” 

/?4 = ^ \-^-B\au + onaal > Ps = 


= 5 ( 7 B. + A (# 2 b 2 + «,b, ) - 


At 


15A 


B , 


A> = 


&2<^2 


where 7 has been introduced, and is given by 

7 = _a* — T— ^ 

' 6>! A V n, A / 
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The terms ao, a x , a 2 and <13 are given by 

A3 2 A3 1 3A3 

a„= T A •, a, = T A J— , a 2 


M 2 
a ej ' 


2 A 3 

a3 = T a \xf x 


a \xe 

For brevity j and a 2 have been introduced, and are given by 

“i 3 = ~ ttj+i'j and oc 2 = (0 • 

Now applying no-slip conditions at the wall it is found that 

i(9 x P 2 +9 2 P x ) = A3£ 2 Ai' + A^BjAi"" +a u B l Ai , {, 

ou\ 

>(«iP 3i +«3iP.) = AJBjtAi; + aI^b,aC' + ;a-» - n 3£ A b, ac 

i(6,P, +« 2 P 2 + ftp.) = A* + A*B,^. 

i=0 ^ 

Matching with the main deck yields 

B 2 k - «2Aio - cvn^iAio = A ^4 2 - P X I - , 

•®3L« + Aig — — “ — ^3L^ B\ Ai 0 = A — Pi I — , 

6 00 

Ba-c-A-l ^ftAi'-'-A* Jfd( = X (a, -Pj-^hp.1- h-M + |a, - 

1—1 (o 

(3.6) 

The integrals occurring in the expressions (3.6) are given by 

i = f_u B 2 (F, ) dFj , and I — J- u B z (F, ) dF, , 

where A is arbitrary and non-zero. The horizontal bar denotes that only the (Hadamard) 
finite part is to be retained. We note that for small r/ the error function can be ap- 
proximated by, A {rj — r/ 3 /3 + T] 5 / 5), and hence / is given by 


J~-l i — 

A 2 1 


85A - 12A 3 


A 26 (30 - 10A 2 -(- 3A 4 ) 


+ 


1 f 43 
26 J 30 - 


435 - 90t/ 2 


IOt ; 2 + 3 r/ 


4 dT 1 ’ 
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with a similar expression for I. In the numerical calculations the quantity H 200 — H] 00 — 
I is found to be —0.29512(2^), this is required in the second order eigenrelation given 
by (3.8 b). 

§3.6 The upper deck solutions and eigenrelations 

The pressure equation in the upper deck yields the bounded solution 

Pi =Pie~ e '\ 

where P\ = P\ (t) is an unknown amplitude function. Now matching with the main 
deck as y — > 0 yields 


P\ — P\ i Q\P\ — A\ (= iOxv^oo ) , 

where v^oo denotes the values of V\ (the main deck normal velocity) as the decks 
transverse coordinate (Y) tends to infinity. This provides the leading order eigenvalue 
problem 

i6\n = AA» Aig. (3.7) 

This relationship determines the wavenumber in terms of the frequency and the skin 
friction A. The higher order pressure functions may be determined, and are given by 

p 2 = (P 2 -0 2 P,y) e-*’*, 

p3L = (Pil - 0 3L Piy ) e~ ei9 , 

Ps = (P 3 ~ ( 03 A + 0 2 P 2 ) y + ^Piy 2 ) e~ d ' 9 . 

To match with the main deck, the unknown functions must satisfy 

p 2 = p 2 -A 1 e 2 l H l0Q , 

.02 

i9\P\H 2oo + 2ifti A\ — i9\A 2 = iA\ 0 2 — iP 2 — i — P\, 

0i 

P3L = P3L , 

-0i Azi — i6nA\ - iPzL - 1 , 

VI 

P3 = P3 + 0 2 Pi H300 — 0i (0i A 2 + 9 2 Ai) I fioo + 2ft j 0i ^ 4 1 P400 > 
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i(e 2 Pi + 0 1 p 2 )h 2oo - ie\A x H 5oo + in 1 P 1 Heoo + *12 1 (a 2 - 

+ *02^4i — iO x A 3 = — ((fij — 0 2 ) v 2oo + (^2 — # 3 ) vioo) — iPz — t - 3 . 1 — a • 

"1 #1 01 

The higher order eigenrelations may now be obtained, 

£>02 _ o 2 0, = *0* (h 2oo - H loo - f) + 2*12! , (3.8a) 

^03L - «3L01 ^ 7 ? = 6^1, (3.86) 

kA5 A 4 

D6 3 - Q. 3 8i ° j = i ( c.t. s + c.t. \ 3 + n.t.) , (3.8c) 

kAs 

where D and D are given by, 



In (3.8) c.t s denotes the terms that arise in conventional analysis; these are exactly 
the same as those occurring in the third order terms of Smith (1979 a), as the non- 
parallelism has no effect at this order. These terms can be directly determined from 
(3.15 b) of Smith (1979 a). The terms c.t • a 3 are those associated with the second 
order basic flow terms and n.t. are the terms that arise due to the t-derivatives (that 
is unsteadiness of the modes). Discussion of the solution of these eigenrelations is 
included in section 5. 

§4 Nonlinear theory for larger disturbances 

We start by setting up the nonlinear problem for an unsteady triple deck. As in 
the linear problem the scalings x = e 3 X and t = e 2 T are introduced. In the Rayleigh 
layer under consideration at the time station t, the flow field in the absence of any 
disturbance is given by 

U = UB(t,Y) + ... , 

where « B (f,K ) is given by the error function. It is also worth re-iterating at this 
point that as Y — + 0, 


u B (<,r)- A(t)T + A 3 (t)F 3 + ---, 
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where A ( t ) = (7rt)~ T . 

The properties of the entire problem may now be determined by examining the 
lower deck, the same transverse coordinate is used as in section 3.2, that is y = e 5 y, 
where y — 0(1). In the lower deck the flow is taken to be of the form, 

u = (eU (X, y, Z, T; t) , e 3 V (X, y, Z, T; t) , eW (X, y, Z, T; t) , e 2 P (X, Z, T; t)) . 

This flow includes both the basic flow and any disturbance quantity. The fact that 
the pressure P is independent of y is given by the wall normal momentum equation. 
The governing equations in the deck are now found to be 


dU 

dT 

dW 

dT 


+ U 


r 

dX 

dW 

dX 



T dU 

dP 

d 2 U 

+ V 

dy 

+ w — = - 

dZ 

dX + 

dy 2 ' 

T r dW 

,rrdW 

dP 

d 2 W 

+ dy 

+ W - 

dZ 

dZ + 

dy 2 ’ 


dU dV dW 
dX + dy + dZ ~ ’ 


(4.1a) 

(4.16) 

(4.1c) 


with boundary conditions 


u = V = W = 0 


at 


y = o, 


(4. Id) 


and, 


U~ \(y + A(X,Z,T)) as 


y 


oo. 


(4.1e) 


Here it is assumed that the spanwise component of the disturbance is confined to the 
lower deck. From the upper deck equations we see that 


V 2 p = 0, with p = P (X, Z, T) and 


dp _ d 2 A 
dy ~ dX 2 


at y = 0 


where y is the transverse upper deck coordinate, and p is a function of X, y, Z and 
T, and is bounded as y — > oo. Now consider a small perturbation to the basic state 
( U — Ay, V = W = P = A = 0) in the lower deck. The perturbation to this system is 
taken to be of order ItCl. In Hall & Smith (1984) initially a bi-modal analysis was 
considered and subsequently a multiple mode analysis, here a single mode is considered. 
It is assumed that h is larger than any positive power of e so the governing equations 
can be assumed to be given by (4.1) for all orders of h. In section 6, the problem 
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involved in making h ~ c" is discussed, where unsteady effects occur at the same order 
as nonlinear effects. The flow is now perturbed and taken to be given by 

U = Xy + hUi + h 2 U 2 + h 3 U 3 + ... , 

V = hVi + h 2 V 2 + h 3 V 3 + ... , 

W = hW] + h 2 W 2 + h 3 W 3 + ... , 

P=hP,+ h 2 P 2 + h 3 P, + . . . , 

A = hA\ -J- h 2 A 2 h 3 A% -|- . . . , 

P = hp\ + h 2 p 2 + h 3 p 3 + . . . , 

The disturbance is considered to be proportional to E and its integer multiples, where 
E is given by 

E = exp (i ($X + 7 Z- SIT)) , 

0, 7 and are real constants for neutral stability, and Z represents the span wise 

_ 3 

variation on the 0(R e 8 ) scale. For any spanwise wavenumber 7 a pair 9 and Q. may 
be determined using the first order eigenrelation given by linear theory. At this point 
a new time coordinate is introduced, namely T defined by 

T = ~T = h~ 2 f. 

X 2 1 

This necessitates the inclusion of further terms in the skin friction in the neighbourhood 
of the neutral time to, as mentioned in Smith (1980); the correction to Smith (1979 
b), so 

t = t 0 + h 2 t, A = A + h?\ 2 + 0(h 3 ), Sl = Sh+h 2 n 2 + 0 (h 3 ), 

where A2 = t(d\ / d€)\t=t 0 . It is necessary to include the further terms in the fl form, 
to satisfy the third order eigenrelation. 

Now substitute the disturbance into the governing equations including the neces- 
sary multiple scales approach for temporal derivatives. The equation of continuity at 
successive powers of h becomes, 

Lx (U i ,Vi,W 1 ) = 0, Lx(U 2 ,V 2 ,W 2 ) = 0, L a (U z ,V z ,W z ) = 0. 

The streamwise momentum equation at successive orders yields, 

L 2 (U\,V\,P\) = 0, 
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l 2 (u 3 ,v 3 ,p 3 ) = -^-\ 2 y^-v,\ 2 -u,^-u 2 

-v^-v^-w^-w™' 


ar ar 

The spanwise momentum equation become, 


az 


az 


at/, 

ax 


l3(Wi,Pi)«0, 


aw-i dWt aty, 
i 3 (w 2 ,p 2 ) = -Ux W- - vi- W- - wk 


ax 


ar 


az ’ 




ar 


ax 


ax 


ax 


1 ay ay az az 

The operators £,, L 2 and £3 are defined by, 

r /rrTrTI „ dU dV dW 

L, (£/, y, W) - Q X + QY + dz ’ 

dU dU dP d 2 U 

L 2 {U,V,P) = — A XY — + XV + — - 

r / xtt 7-1 \ < 9 ^ XT ,avy dp e 2 w 

Lz ( w ' p ^~ ot + XY ax + az ay 2 ' 

The boundary conditions at successive orders are given by 


Ui = Vi = Wi = 0 at 

U\ -» AA, (X, Z, T) 
[7 2 ^AA 2 (X, Z, T) 
U 3 -^XA 3 (X,Z,T) + X 2 A 1 ( 


y = 0 

(* = 1,2,3), 

as 

y -» 00, 

as 

y -*• 00, 

,Z,T) 

as y — » 00 , 


It is worth noting another significant discrepancy between this work and that con- 
sidering a spatially growing boundary layer, that is the slow time derivative does not 
come into the upper deck governing equations, as was noted in the linear section. Now 
consider U\ to be wavelike proportional to E and its conjugate denoted by E. The 
first order disturbance field is now given by 

U\ = tJx (y, t') E + c.c., V\ = V\ ^y, T^j E + c.c., W\ = W\ ^y , Tj E + c.c., 
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P\ - P\ (r'j E c.c., p\ =p n (y,Tj E + c.c., Ai = Ai (f'j E + c.c., 

where c.c. denotes complex conjugate and the disturbance quantities with tildes are 
independent of , Z , and T. By combining the first order equations the two governing 
equations are obtained 




d? 


where 


At OX ’ 

£ = A* (y- 


ff'lh ,9U, A 7 i ti , 

“C-S7- = — rr^i> 


0 A j 




A = *A0. 


At 


(4.2a, 6) 


(4.2c) 


At this point it is advantageous to introduce a further disturbance quantity that is 
Xij ■ The equations (4.2 a) and (4.2 b) may be manipulated and it is found that Xu 
satisfies 


d 3 Xn _ An _ _ 
S u, 


where 


ae 3 




Xu = iOUl + ijWi. 

The solution of (4.3 a) together with a bounded condition at infinity yields, 

XI 1 =B, (r) Ai(£). 

The three dimensional eigenrelation may now be obtained, namely 

A 2 Ai' (£o) = (iOX ) 5 (9 2 + q 2 ) 2 k, 
as given in Hall & Smith (1984), where 


0,1 f 

Co — T and K — I Ai (q)dq. 

(«A)» I 


(4.3a) 

(4-36) 

(4.3c) 

(4.4a) 

(4.46, c) 


The neutral solution of the eigenvalue problem is given by 

£o = -2.298i3. 

In this relationship 9 is real and (4.4 a) yields 

9 » (9 2 +7 2 )* = 1.001A5, 

and, 

Oi = 2.29903A3. 


Thus for a particular spanwise wavenumber 7, a pair 6 and f), maybe obtained such 
that neutrality is assured, a plot of the variation with 7 of 9 and is given in figure 
3. 
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Figure 3 : Variation of 9 and fti with 7. 


The solution of the second order problem is now considered, the disturbance 
quantities now contain second harmonic terms, the disturbance is taken to be of the 
form 


U 2 = EU 2 1 + E- X U$ + E 2 u 2 2 + E-*U&> + 1720 , 

with similar expansions for V 2 ,W 2 , P 2 , p 2 and A 2 . The terms proportional to E 2 are 
considered and it is found that ^1/22, V22, IV22, /*22^ satisfies the system 

dV 22 

2iOU 22 + + 2 ijW 22 = 0 , 

oy 


-2 jt( c ) 


-2iClU 22 + Xy2i6U 22 + AV22 + 2i6P 22 


d 2 u 22 , a ~ 2 ir dih . T ~ xr fT 

— — = —i9U{ - Vi ~ i'iW\U \ , 


- 2 iClW 22 + Xy2i9W 22 + 2 :i 7 P 22 - = -*W, W, - V,^- - iyW 2 , 

along with boundary conditions 

U 22 — V 22 = W22 = 0 at y = 0, 

U 2 2 — > AA22 as y —* 00. 
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Again at this point it is advantageous to re-introduce Xij , this time X22 defined by 


X22 = 2WU22 + 2*7 W22, 


so the system becomes, 

d 2 X22 


r, c d X 22 ( dxu , T/ ^X 

2 £^7" = - X11-ST- + Vi 






ae 


f) 


The solution of which necessitates the introduction of £, given by £ = Note that 

it is not necessary to solve the equation to determine (U21 ? V21 > W21 ? ^21^ as it can be 
seen that the governing equations for these quantities are the homogeneous equations 
( 4 . 2 ); the solution is a linear multiple of the first order sohition. Now finally at this 
order the mean adjustment components are determined by 


dV 2 o _ 
dy U ’ 


d 2 U 20 _ ~ du[ c) fA c)dU, 


XV 2 


20 


dy 2 

d 2 W 20 


= -Vi 


dy 


-v; 


dy 


- *7 (tV, (c) 


Ui - w x u 


; c) ) , 


= —i 0 


dy 2 

with boundary conditions 


(wl 


(c) 


Ui -WiU. 


P) - 




(c) 


dy 


Vi 




V20 = V20 = = 0 


at 


y = 0 . 


By the usual substitution, that is X20 = i&U 2 o + *7^20 ? it is shown that 

dV 2 o 


dt 


= 0 , 


A — — = A’ Vi — - A^ V - Xn + A* — (xiiXn)- 

d £ 3 1 d ( 2 1 d£ 2 j 

Finally the third order equations are considered and again by elementary analysis it is 
known that the third order quantities contain terms proportional to E 3 , E 2 , E 1 , E°, 
E x , E 2 and E 3 , however, determination of the functions associated with E is sufficient 
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to give the required amplitude equation. The functions \ fj31, V31, Wz\,Pz\j are given 
by the solution of the system 

_ d_Un _ _ 7 ^ W ( A -t e M + _ A -*7 

f 3 s * 2 A KK 31 ' * 2 l ^ ^ n/- ' ' ~ ~ ^ 


< 9 ^ 


^ At 3 T 0 
1 d 


K 2 

A3 




(*'* (^ 22 ^ 


(c) 


+ UM 


ae a# ae 


+Ai[ v^ + v^ + vJ^ + v 10 ^ 

' 1 d( d( di di 


+*1 (' 


2 U ' 22 wl c) - W 22 u\ c) + W20 Wi)) , 


d 2 W 31 

a £ 2 


1 aWi *7 ~ A 2 i 0 ( 1 flj A ~ 

_W,= Allr + aT 5 ,+ aH ' ?+ wJ^' 

+ A (2W' 22 & 1 (e) - u 22 w\ c) + u 2 „w,) 

A 3 ^ ' 

J_ f + * ^ + + v 20 Ml 

Ail 1 ^ 1 22 a£ a£ 


+ (w- 22 Wj (c) 

As v 


+ W 20 Wi 


)■ 


with boundary conditions 


U31 = V31 = W31 =0 at y = 0, 


U31 — ► AA 31 + A 2 J 4 n as y -»• 00. 

We make the substitution X31 = i 9 Us\ + 17W31 and it is found that 
a 3 X 3 i 


a ^ 3 




-l ^Xu 

di 


+ I ^ n ^_A*V,^!L-— Vi^^+A*V 3a ^i 

+ 1 2 x a ^ 2 2 1 a ^ 2 a ^ 2 


+Ai^ix« 


A i axn A 


We require a solvability condition for this equation (with the upper deck terms), which 
yields 


dBj 

dT 


— o\B\ (t 2 B\ |£i| 2 , 


( 4 . 5 ) 
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where B\ is the lower deck amplitude function, and where <r 2 is an integral given 
in the appendix (equation (A.l)). Note that o\ may be thought of as a frequency 
perturbation, and if o\ = ift, we may determine ft from the eigenrelation (4.4 a), 
along with (4.4 b,c). If ft is replaced ft + /i 2 ft and A by A + h 2 A (ft = ft 2 and A = A 2 ), 
we obtain the relationship, 

ft aOw-I *oAi + §(;*)* (0 2 + 7 2 )*) 

n = T“ / i i \ -■ (4.6) 

“ A ((i0) 5 (0 2 + 7 2 )’ -A§6>) 

Note that it is not necessary to calculate the value of <x 2 as it has the same value 
as in Hall & Smith (1984), where it is been assumed that A is unity. The velocity 
field has been normalised so that ^-|y =0 = 1. These results are for a single mode, 
whereas Hall & Smith (1984) considered two modes, note that the value of a\ r (which 
corresponds to cr 2r here) is independent of the second spanwise wavenumber thus, 


7 


3.9 

3.6 


m 



2.1 

1.8 

1.5 


@2r 

BEE 

-0.0049 


BB 




-0.0314 

-0.0515 

-0.0919 



Table 1, 7 against <7 2r , as in Hall & Smith (1984). 


Note that ft = q\ {t — and A = A' (to) (t — t n ), hence the relationship (4.6) 
becomes 

A'(i 0 ) (l.46912At) 

A (t 0 ) (2.297205 (02 + t 2) 5 + 5.266128A5 J 

where a superscript prime represents differentiation with respect to t. It is noted that 
A' (to) < 0, thus ftj > 0 after the neutral time t = f n , as is expected from linear theory. 

Note that by considering smaller disturbances unsteady effects may be incurred. 
It is assumed that h is small, but not so small that the leading order effects described 
in this section are changed. In fact the calculations included in this section are appro- 

_ 1 _ i 

priate to 1 >> h >> R e T5 . The case where h ~ R e 715 is now considered, note that it 
is similar to the case discussed in Hall & Smith (1984), for the spatial problem. It is 
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taken that h = HRe^ , with H of order one, and corresponding to a disturbance size 

_ 3 

of 8 — HR e 17 , and it is found that the nonlinear amplitude equation is given by, 

^ = k i (-t 1 +t + fH- 4 )B 1 +a 2 B 1 \B 1 \ 2 . (4.7) 

dT ' ' 

The coefficients in this equations remain the same as in (4.5) so no further numerical 
work is needed, the discussion of the solutions of this equation are included in section 
6 . 


§5 The neutral stability characteristics 

§5.1 Quasi steady flow stability in the triple deck (linear case) 

Here we find the wave characteristics that can be determined from the eigenre- 
lations (3.7 Sz 8). It should be stressed at this point that triple deck theory is only 
applicable for R e >> 1, however we can compute the neutral wavenumber 0„ and 
compare it with the Orr-Sommerfeld predictions. The eigenrelations (3.7,8a & b) are 
unaffected by the boundary layer’s growth, so the condition of neutrality is realised 
by ensuring that f2j, O 2 and are real. 

To ensure that flj and #1 are real we calculate the value of £0 (defined in section 
3.5), required to make the quantity AiQ/(f3/c)(= ci = 1.00065) real. This occurs when 
£0 ~ —2.298^3 - which yields 


3 ^ 


3 ^ 


ft, = XJCocJ = Affij, 0i = Afc* = A* 0] . 


At this point we recall that Ai = but 0i is a constant so we can think of this 

expression as defining the neutral time, 

<„~-0r f - 

7T 

After this time modes will exhibit temporal amplification, and the weakly nonlinear 
analysis described herein will come into effect as the mode attains a finite amplitude. 
It appears that the higher wavenumber modes will start to grow earliest, of course this 
can be seen from figure 1. These results may also be interpreted in terms of the local 
boundary layer Reynolds number, defined by Rf = 2/ v / 7r(fi? e )^ , where 2\fi/y/n is 
the displacement thickness. As noted in section 2, the Orr-Sommerfeld calculations 
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are affected for t = 1/4 so that the value for the displacement thickness is l/2>/F. We 
can express the neutrally stable wavenumber in terms of this Reynolds number as 


0 


n 



7 r 4 



We now proceed to the second order eigenrelations, and note that 



ci 


Ai 20] Aio 
3 + Kile} 


We require that and 02 are real, and thus we find that 


0 2 = 1.55721t“*. 


In order to determine the logarithmic wave characteristics we note that A 3 = — l/(12v7ri^), 
and we find that 

0 3L = -0.56748*“*. 

Subsequent terms in the expansion are effected by the growth of the layer, hence we 
defer their discussion to the next section. We may combine the expressions for 6 \ , 02 
and 03 /, in order to predict the neutral wavenumber 6 n in terms of the local boundary 
layer Reynolds number. We find that 

0 n = 1.0789JV ( l + 3.37079 Rj) - 5.70354 Rj) In R s * +■••), 

this value is compared with the neutral value found from Orr-Sommerfeld theory, the 
results are displayed in figure 4. In order to predict the upper branch we would need 
to consider a more complicated structure in which the critical layer is distinct from 
the viscous wall layer. 
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Figure 4: Comparison of neutral wavelength, Orr-Sommerfeld and triple deck 

The agreement is reasonable, but it should be noted that analysis relies on different 
powers of the Reynolds numbers being distinct, however at the values of Rg» displayed 

_ I _ 1 

in figure 4, In and R 6 , 4 are not totally distinct. However, the triple deck 
theory does appear to capture the lower branches characteristics. 

§5.2 Unsteady effects on triple deck stability 

In this section we discuss the effects that the boundary layer’s growth has on 
the stability of the Tollmien-Schlichting waves. As has been mentioned several times 
earlier, the evolution of the layer changes the characteristics at one order higher than 
in the spatial case of Smith (1979 a), thus it is perhaps even more crucial in these 
temporal cases that the layers evolution is included in the analysis. This is to be 
expected, since during one spatial wavelength of the wave the boundary does not 
evolve (spatially) as much as it would during one period (temporally). The leading 
order effect of the boundary layer’s growth occurs in the eigenrelation ( 3.8 c), which 
determines O 3 and $ 3 . 

We consider the total temporal derivative of a flow quantity $, and define its 
growth rate as, 

1 d$ 

a ~ $ d v 
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We require that the quantity a is purely imaginary, which for a problem which does not 
exhibit any temporal changes translates to ensuring that the f}’s are real. However, 
we know that d\ —* d t + e~ 2 &r, so that 


a = 


-2 


— ifli — eiQ , 2 — e 2 In tSlzi, — e 




Hence in order for a to remain totally imaginary (to the order considered), we require 

Images) = —Real 

Using this condition together with eigenrelation (3.8 c), we can determine 0 3 and 
[Real (ft 3 ) - Imag (I 1 -)] . 

It is not clear as to which quantity should be used in order to represent the 
disturbance’s growth. In Smith (1979 a) the wall pressure was used, and the results 
which were obtained allowed enhanced agreement with the experiments of Ross et 
al (1970). Smith noted that this implied that the boundary layer’s growth had a 
stabilising effect, however if other quantities were chosen to represent the disturbance, 
the conclusions could be the opposite. In Gaster (1974) and Eagles k Weissman 
(1975), it was noted that the difficulty is in defining the instability in some absolute 
sense, however this method should allow good agreement with a particular experiment 
(in which the disturbance is measured using a certain quantity). We are not aware of 
any experimental data concerning the susceptibility of a Rayleigh layer to Tollmien- 
Schlichting waves. If, for example we take the wall layer normal velocity (Vi ) to be 
representative of the disturbance, then the layers growth is seen to have a destablising 
effect. 

The reason for the enhanced agreement with experiment seen in Smith (1979 a), 
is likely to have its origin in the size of the correction, it is now of order R e ® which 
is larger than that of Bouthier (1973) and Gaster’s (1974) (order R e 5 ) successive 
approximation technique. As mentioned previously our correction is of order R e * 
which represents an even larger change, so we expect even better agreement between 
the theory and experimental results. 



27 



§6 Discussion of the weakly nonlinear stability. 

The relationship between the equilibrium amplitude solutions of (4.5) and those 

of the unsteady flow amplitude solutions of (4.7) is now considered. The latter as 

!_ 

previously stated corresponds to a disturbance of size h ~ R e 16 . The former is for 
larger disturbances where h > R e . Thus we are required to solve, 

^- = iClB + cj 2 B\B\\ (6.1) 

at 

and, 

~ = iqi (-<» + H~ a t) B + <j 2 B\B\\ (6.2) 

where B = B\ is the lower deck amplitude function as in (4.3 c). The first equation 
(6.1) has the stable equilibrium solution, 

I B| 2 = (6.3) 

0 2r 

which bifurcates supercritically at fb = 0 from the zero solution for increasing Reynolds 
number, where <r 2r is always less than zero as shown in table 1. Increasing the Reynolds 
number corresponds to progressing through time, thus increasing to- The zero solution 
is unstable beyond Cl equals zero, so any perturbation to this tends to the solution (6.3). 
Now note that sufficiently close to the neutral time t n (specifically within h ~ R e ), 
the appropriate equation is (6.2). Drawing on the conclusions of Hall Sz Smith (1984) 
we let B = p exp where p and 4> are real functions of T. It was shown that the 
large T behaviour of p is given by, 

2 q i H~ 4 T f qp ^ 

<7 2r k / 

It may be seen from the numerical results of section 5 that q, and <r 2r are both always 
negative. This equation implies that the large T behaviour may be inferred from (6.2) 
neglecting the dB / dt term. This equilibrium solution corresponds to letting Cl —* 0 in 
equation (6.1), which corresponds to the disturbance size being decreased. Thus the 
interval over which unsteady terms induce a contribution to the amplitude equation 
is crucial. In this interval a small disturbance ( h ~ Rj 15 ), develops smoothly into a 
finite amplitude Tollmien-Schlichting wave. Only ‘later’ (that is for relatively large 
time) does the amplitude tend towards the amplitude predicted by (6.1). 
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§7 Conclusions 

We have shown that in theory a Rayleigh layer is unstable to Tollmien-Schlichting 
waves. An Orr-Sommerfeld calculation is performed and it is found that beyond a 
certain Reynolds number, disturbances with fixed spatial wavenumber comparable to 
the displacement thickness are temporally unstable. As the Reynolds number increases 
further the structure of the neutral modes can be described by a triple deck method, 
and eigenrelations are calculated and used to determine the wave characteristics. A 
multiple scales technique is used to determine the effect of the growth of the layer on 
the waves’ stabilities. In summary it is found that this depends on which quantity 
is chosen to represent the waves intensity. A Stuart Watson amplitude equation was 
derived in the neighbourhood of the neutral time, and the waves were found to bifurcate 
supercritically from the zero solution. An amplitude is found that represents a crucial 
stage in the development of the wave, this pahse allows the wave to evolve smoothly 
from its linear stage to its equilibrium amplitude. 

Although this work in this paper represents fairly minor modifications to various 
earlier papers; the methods are applied to an entirely different physical problem. It 
will be of interest to determine whether the upper branch modes can also be compared 
with the Orr-Sommerfeld calculation. The calculation will be solved using the methods 
discussed in Bodonyi & Smith (1981), wherein the effect of non-parallelism on upper 
branch Tollmien-Schlichting waves in a Blasius boundary layer is discussed. The upper 
branch problem is far more complex, since the critical layer is distinct from the wall 
layer. Within the context of this spatially homogeneous problem it will be interesting 
to find what acoustic radiation may be generated by these waves, and whether the 
waves on the upper branch or the lower branch of the neutral curve produce the most 
sound. The noise can be calculated using methods similar to those discussed in Tam 
& Morris (1980). 
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Appendix 

The double integral a 2 involved in the weakly nonlinear calculation occurring in equa- 
tion (4.5), is given by 
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a 2 = 


0(* a + 7 2 ) , Aitfo)-Ai"(£ o ) 

where the function G(£) is given by, 
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and the form of F(£) included in the integral <7( a ) is given by, 
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the function H(£) is given by, 


a*ff( Q 

ap 
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where the operators L and L 2 are defined as, 
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The terms c.t.\ 3 occurring in equation (3.8), that is those due to the second order 
basic flow terms, are given by 


c.t.\ 3 = 



£o 


where /(£) as given in section 3.5. 

The terms due to non-conventional analysis in equation (3.8), axe 
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